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Abstract. We discuss a relation between the QCD thermodynamics obtained 
from a statistical analysis of particle production in heavy ion collisions at SPS and 
RHIC energies and recent LGT results at finite chemical potential. We show that 
basic thermodynamic properties obtained from the phenomenological statistical 
operator of a hadron resonance gas that describes particle yields in heavy ion 
collisions are consistent with recent LGT results. We argue that for T < T c the 
equation of state derived from Monte-Carlo simulations of two quark-flavor QCD 
at finite chemical potential can be well described by a hadron resonance gas when 
using the same set of approximations as used in LGT calculations. We examine 
the influence of a finite quark mass on the position of the deconfinement transition 
in temperature and chemical potential plane. 



1. Introduction 

The detailed analysis of particle production in heavy ion collisions has shown that 
in a broad energy range from SIS through AGS, SPS up to RHIC particle yields 
resemble that of chemical equilibrium population 1 . At SPS and RHIC the freeze- 
out parameters, the temperature T/ and the energy density e/, predicted by the 
presently used thermal model of hadron resonance gas, agree well with recent results 
from the lattice on critical conditions required for deconfinement E] ■ The above 
quantitative agreement of freeze-out and critical parameters suggests that at SPS and 
RHIC chemical freeze-out appears in the near vicinity or at the phase boundary 3 . 
If this is indeed the case then the phenomenological statistical operator of hadron 
resonance gas Zhg should also provide a consistent with LGT description of QCD 
thermodynamics in the confined, hadronic phase [2]. Here we show that the basic 
qualitative properties of Zhg resulting from its dependence on T and baryon chemical 
potential fiB are also present in recent lattice results of two flavor QCD at finite 
chemical potential. In addition, when imposing the fixed energy density condition 
for deconfinement the hadron resonance gas partition function Zhg describes [2] the 
quark mass dependence of the lattice critical temperature at fj,B = and the position 
of the phase boundary in the (T, /ZB) _ plane at small [Ib- 
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Figure 1. Left-hand figure: temperature dependence of AP/(\ q T' 2 ) ratio for 
two different values of fi q /T. The 2— flavor lattice results are from |3J. The lines are 
the hadron resonance gas model values with (broken-line) and without (full-line) 
Taylor expansion of cosh(/ig/T)). The right-hand side shows the temperature 
dependence of the ratios of second, forth and sixth order coefficients in the Taylor 
expansion of thermodynamic pressure. The points are lattice results from 
The lines at T < To are the values obtained from the expansion of cosh(/jg /T)) 



2. Resonances essential degrees of freedom near deconfinement 

The phcnomcnological partition function used in the description of particles 
production in heavy ion collisions was, following Hagedorn, constructed as a non- 
interacting hadronic gas which is composed of all hadrons and resonances. In the 
Boltzmann approximation, suitable for the moderate values of \i B < nijv and T > 50 
MeV, there is a factorization of T and hb/T dependence in relevant observables 
characterizing baryonic sector of the system. The basic quantity is the pressure 
AP = P(T,n B )-P(T,n B =0) 

_ = F(T)(cosh(^ - 1) with F(T) ~ J dmp(m){ — ) 2 K 2 ( — ) (1) 

from which n B and the baryon number susceptibility A\ B are obtained as the first 
and second order derivatives with respect to /is, respectively. 

The obvious consequences of the factorization in Eq.(l) is that any ratio of n B , 
AP and A\ B at fixed fJ- B /T should be independent of T. This is the property which 
can be directly checked with recent LGT results. In Fig.(l) we show as an example 
the ratio of AP/ A(x B T 2 ) for two different values of (J- q /T as function of T* It is 
clear from Fig.(l) that the factorization predicted by Zhg is also seen in LGT results 
for T < To, that is in the confined phase of QCD. 

The second transparent feature of Zhg is that the n B /T dependence appears 
through a cosh-function. The lattice results at finite fi q were obtained using a Taylor 
expansion of the pressure with respect to the quark chemical potential. So far, the 
Monte Carlo results are available for the second c 2 , fourth C4 and sixth cq order 
coefficient in the Taylor series of P with respect to [i q jT . If the lattice thermodynamics 
was consistent with the prediction of Z B q then one would expect C2/C4 = 3/4 and 

* The quark chemical potential in Fig.(l), fi q = t^Hb 
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Figure 2. Left-hand figure: temperature dependence of baryonic pressure 
for different values of fi q /T. The lattice results are from [I]. The lines are 
the statistical model results 0. The right-hand figure shows lattice results on 
phase boundary curve (line with errors) together with phenomenological freeze- 
out values of T and fig (points) obtained from the analysis of particle production 
in heavy ion collisions 0. Short-dashed and dashed-dotted lines are the statistical 
model results obtained under the condition of fixed e ~ 0.6 GeV/fm 3 with 
m v ~ 0.77 GeV and m, ~ 0.14 GeV respectively. Also shown (full— line) is 
the phenomenological freeze-out curve of fixed energy/particle~ lGeV from [?]. 



Cq/c4 = 0.3 being the coefficients in the expansion of cosh(/xs/T) in Eq.(l). Fig.(l- 
right) shows the corresponding ratios obtained on the lattice and the results obtained 
from Eq.(l). Within statistical error the agreement is indeed seen to be justified. 

The pressure calculated on the lattice (see Fig.(2-left) ) increases abruptly when 
approaching deconfinement transition from the hadronic side. If the phenomenological 
statistical operator Zhg is of physical significance then this increase could be due 
to resonance formation. To check the importance of resonances near deconfinement 
one would need to reproduce lattice results on the T-dependence of P at fixed 
Hq/T. However, to quantify this dependence one needs to implement the same set 
of approximations in Eq.(l) as those being used on the lattice. First of all current 
lattice results in Fig.(2-left) are obtained with quite a large quark mass corresponding 
to rn,,. ~ 770MeV. This also distorts the baryon mass spectrum. Its pion mass 
dependence can be deduced from lattice calculations at zero temperature. We use 
the following ansatz for the parametrization of the dependence of baryon masses on 
the pion mass 

m vn 

where A — 0.9 ±0.1, m* is the distorted hadron mass at fixed m n and m is its 
corresponding physical value. 

The lattice results were obtained in 2-flavor QCD, thus there is no contribution of 
strange baryons in Eq.(l). Due to the Taylor expansion of P in the lattice calculations 
one also needs to perform a similar approximation in Eq.(l). In Fig.(2-left) the 
lattice results are compared with Eq.(l). The T dependence of QCD thermodynamics 
obtained on the lattice is seen in Fig.(2-left) to be consistent with the predictions of 
Zhg- 
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The value of T c was shown on the lattice to be dependent on the pion mass 
UJ. The remarkable feature of this dependence is that e at T c is almost constant, 
independent of m ff . This would suggest that deconfmement is density driven and can 
be obtained from the condition of fixed energy density. Fig.(2-right) shows recent 
lattice results on the position of the phase boundary line in 2-flavor QCD within 
the Taylor approximation and with the quark mass such that m n ~ 770MeV. The 
condition of fixed energy density e ~ 0.6 GeV//m 3 with e obtained from Eq.(l) is 
seen to coincide with lattice results. Decreasing the pion mass to its physical value 
and including complete set of resonances expected in (2+l)-flavor QCD results in 
the shift of the position of the phase boundary line towards phenomenological freeze- 
out condition of fixed energy/particle~ lGeV. The splitting of freeze-out and phase 
boundary line appears when the ratio of meson/baryon multiplicities reaches the unity. 

3. Conclusions 

We have shown that the statistical operator of hadron resonance gas used to describe 
particle yields in heavy ion collisions provide satisfactory description of recent lattice 
results on QCD thermodynamics at finite chemical potential. In particular the basic 
property of this operator like e.g. factorization of temperature and chemical potential 
dependence is obviously confirmed by the lattice results. The ratios of the coefficients 
in the Taylor expansion of thermodynamic pressure are also well described by the 
expansion of cosh-function predicted by the hadron resonance gas. These results are 
independent from the particular choice of the quark mass in the lattice calculations 
and are also to large extent free from lattice artifacts. The phenomenological partition 
function was also shown to describe quantitatively the temperature and chemical 
potential dependence of basic thermodynamical observables bellow deconfmement. 
This indicates that hadron resonance gas partition function is a good approximation 
of QCD thermodynamics of the hadronic phase. That is why this partition function 
is also successful in heavy ion phenomenology. 
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